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ABSTRACT 

We study the optical properties of the host galaxies of nuclear 22 GHz (A = 1.35 cm) water masers. 
To do so, we cross-match the galaxy sample surveyed for water maser emission (123 detections and 
3806 non-detections) with the SDSS low-redshift galaxy sample (z < 0.05). Out of 1636 galaxies with 
SDSS photometry, we identify 48 detections; out of the 1063 galaxies that also have SDSS spectroscopy, 
we identify 33 detections. We find that maser detection rate is higher at higher optical luminosity 
{Mb), larger velocity dispersion (a), and higher [O III] A5007 luminosity, with [O III] A5007 being the 
dominant factor. These detection rates are essentially the result of the correlations of isotropic maser 
luminosity with all three of these variables. These correlations are natural if maser strength increases 
with central black hole mass and the level of AGN activity. We also find that the detection rate 
is higher in galaxies with higher extinction. Based on these results, we propose that maser surveys 
seeking to efficiently find masers should rank AGN targets by extinction-corrected [O III] A5007 flux 
when available. This prioritization would improve maser detection efficiency, from an overall ~ 3% 
without pre-selection to ~ 16% for the strongest intrinsic [O III] A5007 emitters, by a factor of ~ 5. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — radio lines: galaxies — 
masers 



1. INTRODUCTION 

Water maser emission at 22 GHz (A = 1.35 cm) 
is currently the only tracer of warm dense molec- 
ular gas in the inner parsec of active galaxy nu- 
clei (AGNs) and has been detected to date in more 
than 100 AGNs (e.g., " 
2005 ^ iKondratko et al.1 



Braa tz et al.lll996t iHen kcl ct al. 
2006U iBraatz fc Gugliuccil T2008: 



Greenhill et al. I I2008IL Some of these masers are as- 



sociated with rotating, highly inclined disk structures 
close to the central engines ("disk masers") and have 
been used for a broad variety of astrophysical stud- 
ies, including the mass estimation of supermassive black 
holes, the mapping of accretion disks, and the de- 
term i nation of geometric di s tances (e.g., iMivoshi et al.1 
19951: [Greenhill et al.l fl997L IGreenhill fc Owinnl Il997 



ferred X-ray luminosities are subject to large uncertain- 
ties owing to high columns, and the true (i.e., beamed) 
maser luminosities are unknown in most cases. 

The overall detection rate of nuclear water masers is 
only ~ 3%. Even if AGNs with higher X-ray luminosity 
and/or higher obscuring column more likely host masers, 
there is no existing large sample of AGNs with X-ray 
data available for target selection. However, if masers 
are preferentially found in galaxies with certain optical 
properties, we can improve maser detection efficiency by 
selecting galaxies with these properties as targets from 
existing large galax y surveys, such as the Sloan Digital 
Sky Survey (SDSS; I York et all l2000h. the 2d F Galaxy 



Redshift Survey (2dFGRS: IColless et al 



6dF Galaxy Survey (6dFGsT Tjones et aT 



2004) 



2001), and the 



Ishihara et all 120011: IGreenhill et alJ 12003b IBraatz et al] 
20101 : IKuq et al.ll2010f >.~ 

Nuclear water masers have been claimed to be associ- 
ated with Seyfert 2 or low-ionizatio n nuclear em i ssion- 
line region (LINER) sy stems (e.g.. IBraatz et all 11997b 
IKondratko et alj |2006b). It is also plausible that AGNs 
which host masers are more likely associated with high 
X-ray obscuring columns (Nh) than those without maser 
detections (e . g., IBraatz et al.lll997tlMadeiski et al.ll2006t 
Zhang et all 120061: IGreenhill et al.l 120081: iZhang eTaLl 
20101) . There also appear to be correlations of 
isotropic maser luminos ity with the X-ray luminosity 
(IKondratko et all l2006bH an d the far-infrared (FIR) lu- 
minositv (iHenkel et al.ll2005h of the host AGNs, though 
the underlying mechanisms are not clear. As additional 
words of caution in interpreting these correlations, the in- 
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The goal of this work is to systematically investi- 
gate the optical properties of maser host galaxies. We 
cross-match the SDSS low-redshift galaxy catalog with 
the complete galaxy sample surveyed for maser emis- 
sion. We find that maser detection rate is higher 
at higher optical luminosity, larger velocity dispersion, 
higher [O III] A5007 luminosity, and higher extinction. 
We present these results in Section [2l In Section [3j 
we suggest that a plausible explanation of these results 
is that maser strength is correlated with the central 
black hole mass and the AGN activity of the host galax- 
ies. In Section 2J we suggest that maser surveys rank 
AGN targets by extinction-corrected [O III] A5007 flux, 
which should greatly improve the detection efficiency. 
We adopt a ACDM cosmology with Q m = 0.3, fi A = 0.7 
and H = 70 km s _1 Mpc" 1 . 



2. DETECTION EFFICIENCY 
2.1. Data 

2.1.1. The complete galaxy sample surveyed for maser 
emission 
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Fig. 1. — Top panel: Color-magnitude diagram. For comparison, 
we show the distribution of the whole low-z photometric sample as 
the gray scale. The contours enclose 40%, 80%, and 90% of the 
sample. Blue triangles represent non-detections and red squares 
indicate maser detections. We also show the disk masers with 
green diamonds but we cannot draw robust conclusions because 
of the small sample size. Middle panel: Distribution of Mb- We 
show the ratio of the number of galaxies per 0.5 mag to the size of 
each sample. Bottom panel: Detection rate as a function of Mg. 
To guide the eye, we show two horizontal dotted lines at 5% and 
10%. The error bars represent Poisson errors. The detection rate 
appears to be higher at higher luminosity. 

To construct a complete sample of galaxies surveyed 
to date for maser emission, we combine the catalogs 
(as of December 1, 2010) maintained on the website of 
the Megamaser Cosmology Project (MCF0) and that of 
the Hubble Constant Maser Experiment (HoME0). For 
maser detections, we use the MCP catalog that is com- 
plete. We however exclude those masers known to be as- 
sociated with star-forming regions (IC 10, M 33, IC 342, 
M 82, NGC 253, NGC 3359, NGC 3556, NGC 2146, He 
2-10, NGC 4038/39, NGC 4214, NGC 5253), as noted 
in the either HoME or MCP catalogs. Although it is 
not labeled in either catalog, we also exclude NGC 4194, 
since it is an ongoing m erger and has strong star forma- 
tion in the center (e.g., lBalzanol l.198 31) and the detected 
maser may not be associated with its nucleus. There 
are 123 detections in total and we list them in Table [TJ 
Among these, at least 41 are probably associated with 
disk structures, as noted in the HoME catalog. The evi- 
dence is either from direct mapping of the emission distri- 
butions using Very Long Baseline Interferometry (VLBI) 

km s - 

4 https : //saf e .nrao . edu/wiki/bin/view/Main/MegamaserCosmology Project 

5 https : //www. cf a. harvard. edu/~lincoln/demo/HoME/index. html 6 http : //www . mpa-gar ching . mpg . de/SDSS/DR7/ 



Fig. 2. — Similar to Figure [T] but with velocity dispersion (<r). 
Note a lower than 70 km s — 1 is not reliable due to the instrumental 
resolution of the SDSS spectrograph; we do not consider galaxies 
below that limit. The detection rate is higher at larger a. 

or inferred from spect roscopy (e.g. JMadejski et al.1120061 
IGreenhill et al.l I2008D . For galaxies without successful 
maser detections (non-detections), we combine both cat- 
alogs from HoME and MCP to build the whole sample. 
After removing duplicates between the catalogs, we have 
3806 non-detections in total. The overall detection rate 
is therefore about 3%. 

We note that this complete sample consists of galaxies 
surveyed using different telescopes with different detec- 
tion sensitivities. The highest sensitivity co mes from the 
Gree n Bank Telescope (GBT) survey (e.g. JBraatz et all 
12004 . The lcr rms sensitivity of the GBT survey is ~ 3 
mJy per 24.4 kHz (~ 0.33 km s _1 ) channel (jBraatz et al.l 
2004). Assuming a characteristic maser linewidth of 
10 km s~ , this corresponds to a 3cr maser flux limit of 



0.1 Jy km s 

2.1.2. The SDSS low-z galaxy sample 

To systematically study maser detection efficiency, we 
require a complete parent sample. The SDSS survey 
has provided such a sample of galaxies with uniform 
imaging and spectroscopy. For spectral properties, we 
use the measuremen ts by the MPA-JHU groupQ (e.g., 
iTremonti et"aT][200l . We use the latest version that cor- 
responds to SDSS Data Release 7 (DR7. lAbazaiian et al.l 
2009). We choose to look at velocity dispersion (a, in 
) and [O III] A5007 luminosity (£[oiii]A5007i m 
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Fig. 3. — The left panels are similar to the lower two panels in Figure[T]and[2] but with observed [O III] A5007 luminosity (£[oin]A5007,obs 
In the right panels, we correct the extinction to obtain the intrinsic [O III] A5007 luminosity (£[om]A5007,cor) using £[oni]A5007,cor : 
^[Oin]A5007,obs/((Ha/H/3)/(Ha/H/3)o) 2 ' 94 , where Hct/H/3 is the observed Balmer decrement and we assume the intrinsic Balmer decremei 
(Hc?/H/3)o = 3. The detection rate is higher at higher [O III] A5007 luminosity. The effect is stronger for £[om]A5007,cori implying th; 
detections on average have higher extinction than non-detections. 



erg s x ), since velocity dispe rsion is closely related to 
central black hole ma ss (e.g.. iFerrarese fc MerrittJ 120001 : 
IGebhardt etall 120000 and [O III] A500 7 luminosity is 
well-c orrelated with AGN activity (e.g.. iHeckman et al.l 
I2005D . It is well-known that [O III] A5007 can be 
severely obscured by mater ial in the host galaxy (e.g., 
IDiamond-Stanic et al.| [2009V We therefore calculate the 
intrinsic [O III] A5007 luminosity (£[om]A5007,cor) by 
correcting the ob served -kfoiniA5007,ob s using the follow- 
ing formula (e.g., iBassani et al~lll999l ): £[om]A5007,cor = 
i[oiii]A5007,ob s /((Ha/H^)/(Ha/H^) ) 2 - 94 , where Ha/H/3 
is the observed Balmer decrement and we assume the in- 
trinsic Balmer decrement (Ha/H/3)o = 3. Because of the 
instrumental dispersion of the SDSS spectrograph, ve- 
locity dispersion measurements smaller than 70 km s~ 
are not reliable^; we thus only consider galaxies with 
g > 70 km s -1 . Finally, since only four (4C +05.19, 
SDSS J0804+3607, Mrk 34, and 3C 403) out of the 123 
maser detections are farther than z = 0.05 and all of 
them are not in the MPA-JHU catalog, we limit the sam- 
ple to low-redshift galaxies with z < 0.05. At the faint 
end, the flux limit of the SDSS spectroscopic survey is 
r = 17 .77, which corresponds to Mb ~ —18 at z = 0.05. 

Due to the difficulty of automatic photometric process- 
ing of big galaxies, the SDSS catalog is missing many 

7 http: //www. sdss . org/dr7/algorithms/veldisp .html 



nearby, bright galaxies, even though they are contained 
within the SDSS imaging footprint. For photometry, 
we therefore use the low-z catalog (z < 0.05) from 
the NYU Value Ad ded Galaxy Catalog (NYU-VAGC0; 
IBlanton et al.|[2005l) . This low-z photometric catalog in- 
cludes any low-redshift galaxies fro m the Third Reference 
Cata l og of Bright Galax ies (RC3; Ide Vaucouleurs et al.l 
[19911 ICorwin et al1ll99l for which we have ugriz imag- 
ing from SDSS, but which are not in the SDSS catalog. 
We use the latest version of th is catalog that corresponds 
to SDSS Data Release 6 fDR6. IAdelman-McCarthv et al.l 
120061) . We have compared the photometry of those galax- 
ies in the SDSS catalog with that from DR7 and found 
they are very consistent, therefore using DR6 for pho- 
tometry should not introduce any bias. We derive ab- 
solute magnitudes using the kcorrect package (v4.1.4; 
IBlanton fc Roweisl [2007T ). For easier comparison with 
previous studies, we choose the B band magnitude Mb 
to indicate optical luminosity. Note here the magnitude 
is the total magnitude for the whole galaxy. 

We cross-match maser detections and non-detections 
with the SDSS low-z sample, and identify 48 detections 
(15 disk masers) and 1588 non-detections with SDSS 
photometry, among which, 33 detections (10 disk masers) 
and 1030 non-detections have reliable spectral measure- 

8 http: //sdss .physics .nyu. edu/vagc/ 
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Fig. 4. — Similar to the lower two panels in Figure [T] and [2] but 
with Balmer decrement Hct/H/3. The detection rate is higher at 
higher Ha/H/3. 




log 10 ([NO] X6584/HJ 

Fig. 5. — Emission line diagnostic diagram (Baldwin ct al. 1981, 
BPT). For comparison, we show the distribution of the whole low- 
z spectroscopic sample as the gray scale. The contours enclose 
40%, 80%, and 90% of the sample. The dotted and dashed lines 
are the demarca tion lines separa t ing A GN a nd star-formin g galax- 
ies defined by iKauffmann ct al. (2003) and Kcwlcy ct al. (2001). 
The solid vertical and horizontal lines at [N II] A6584/Ho=0.6 
and [O III] A5007/H/3=3.0 are the conventional separating lines 
for Seyfert 2 galaxies (above) and LINERs (below). There are 25 
(8) detections above (below) the horizontal line, compared to 296 
(734) non-detections. 



ments, i.e., with ct > 70 km s~ . We note that the over- 
all detection rate in this sample is ~ 3%, the same as in 
the total sample, so using this sample should not intro- 
duce a bias in the analysis presented below. We present 
measurements for detections in Table [1] In the next sub- 
section, we study the dependence of maser detection effi- 
ciency on optical luminosity and spectral properties using 
these samples. 

2.2. Results 
2.2.1. All galaxies surveyed for water masers 

Figure [TJ O and present maser detection efficiency 
as a function of Mb, a (in log scale), £[om]A5007,obs (i n 
log scale), and £[oin]A5007,cor ( m l°g scale). In the top 
panels of Figure [T] and [H we show the detections (blue 
triangle) and non-detections (red square) in the color- 
magnitude/cr diagram. For comparison, we also show 
the distribution of the whole SDSS low-z sample in gray 
scales and contours. Masers are concentrated in systems 
with larger B — R and total luminosity, but there are very 
few blue, low luminosity (presumably disk-dominated) 
galaxies that have been surveyed for emission. We also 
show the disk masers with green diamonds, but we can- 
not draw robust conclusions because of the small sam- 
ple size. In the middle panels of Figure [1] and [5] and in 
the top panels of Figure [3l we show the distribution of 
detections and non-detections. Performing Kolmogorov- 
Smirnov test yields P-values 0.01, 0.02, 2 x 10 -8 , and 
3 x 10~ 9 , respectively, indicating a low probability that 
they are drawn from the same distribution, particularly 
for i[oni]A5007- The bottom panels present the detection 
rate, which is apparently higher at brighter Mb, larger tr, 
and higher £[om]A5007- The dependence of the detection 
efficiency on a appears stronger than that on Mb, while 
the dependence on £[oin]A5007 is more striking than that 
on both Mb and a. 

Figure [3] also shows that maser detection efficiency 
depends more strongly on £[om]A5007.cor than on 
-k[Oin]A5007,obs; implying that detections could on aver- 
age have higher extinction than non-detections. Previ- 
ous studies show that AGNs which host masers are more 
likely associated with high X-ray obscuri ng columns than 
those without maser de tections (e.g., IGreenhill et al.l 
I2008t iZhang et al.ll2010T) . We therefore investigate the 
detection efficiency as a function of the observed Balmer 
decrement in Figured] In the top panel, a Kolmogorov- 
Smirnov test yields a P-value 0.009, indicating that de- 
tections and non-detections are likely drawn from differ- 
ent distributions. The lower panel shows the detection 
rate is indeed higher at higher extinction. 

In Figure El we show the em ission line diagnos- 
tic diagram ([Baldwin et al.l 119811 BPT). We show 
two widely-used demarcation criteria that separate 
AGNs (t o the right) and sta r-forming galaxies (to the 
left), by iKewlev et all (]2001l the dashed line) and by 
IKauffmann et all (120031 the dotted line). The vertical 
solid line at [N II] A6584/Ha= 0.6 and the horizontal 
solid line at [O III] A5007/H/3= 3.0 are the conventional 
demarcation lines for Seyfert 2 (above) and LINE R-likc 
(below) galaxies (e.g.. lVeilleux fc Ost erbrock 1987). It is 
apparent that most surveys have mainly targeted AGNs. 

Among AGNs, masers are clearly more often detected 
in Seyfert 2 galaxies rather than LINERs. For exam- 
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Fig. 6. — The same as the lower two panels in Figure [T] but 
for galaxies with £[oin]A5Q07,cor > 10 40 ' 5 erg s — 1 only. Note we 
double the binsize because of a smaller sample size. Among the 33 
detections, 25 have £[oin]A5007,cor > 10 40 5 erg s — 1 ; and among 
1030 non-detections. 242 have £[om]A5007,cor > 10 40 5 erg s _1 . 
The overall detection rate (~ 9%) is higher than that (~ 3%) 
without prc-sclection in Figurefl] The detection rate is also higher 
at higher luminosity. 

pie, in Figure[5j there are 25 maser detections associated 
with Seyfert 2 galaxies (out of 321 in total, above the 
horizontal line), but only eight associated with LINERs 
(out of 742 in total, below the horizontal line). It is 
yet unkown w h at fraction of LI NERs are AGNs (e.g., 
IHo et all [2001 ISarzi et all l2010h and we are not sure 
about the underlying mechanisms responsible for the 
preferred association of masers with Seyfert 2 galaxies 
over LINERs. However, if LINERs are low-luminosity 
count erparts of Seyf ert 2 galaxies, as some studies claim 
(e.g.. IHo et al. 2003]), then the low maser detection effi- 
ciency among LINERs could be simply a reflection of the 
L[Oin]A5007 dependence we found in Figure [31 since the 
majority of LINERs without maser detections (~ 90%) 
have i[oni]A5007,cor < 10 40 - 5 erg s _1 . After all, LINERs 
are known to have low [O III] A5007 luminosities relative 
to Seyfert 2 galaxies (e.g.. iHeckman et"al~ll2004[ ). 



2.2.2. Galaxies with £[oin]A5007,cor > 10 erg s 1 only 

We have shown that maser detection rate depends 
more strongly on %pni]A5007 than on Mb and a, as in 
Figures [TJ [2j and [3] However, these three quantities 
are themselves correlated in the galaxy population. We 
here investigate further whether there is a residual de- 
pendence of detection efficiency on Mb and a even for 
galaxies with strong [O III] A5007 emission. 
Among the 33 detections, only eight have 



out of 

10 40.5 



[OIII]A5007,cor 

' 1030 



< 10 40 ' 5 ergs 
non-dctcctions have 



erg s 



Meanwhile, 789 

i[OIII]A5007,cor < 

We therefore pre-select galaxies with 



Fig. 7. — The same as the lower two panels in Figure [2] but for 
galaxies with £[oin] A5007,cor > 10 40 ' 5 erg s — 1 only. The overall 
detection rate is higher than without pre-selection in Figure [2] 
The detection rate is higher at larger a. 
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Fig. 8. — The same as the lower two panels in Figure [4] but for 
galaxies with £[oin]A5007,cor > 10 40 ' 5 erg s — 1 only. The overall 
detection rate is higher than without pre-selection in Figure [4] 
The detection rate is higher at higher Ha/H/3. 
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242/1030 non-detections) and show the detection rate 
as a function of Mb and a for this pre-selected sample 
in Figure [6] and [7] We also show the detection rate as 
a function of Balmer decrement Ha/H/3 in Figure 
Compared to the results for the whole sample (Figure 
[U [2 and SI), there are still detections over the whole 
range of Mb, cr, and Ha/H/3, and the detection rate 
is still higher at brighter Mb, larger a, and higher 
Ha/H/3. The overall detection rate, however, is ~ 9% 
compared to ~ 3% without pre-selection. We therefore 
conclude that among Mb, cr, [O III] A5007, and Ha/H/3, 
[O III] A5007 is the dominant factor when determining 
maser detection efficiency. 

3. DISCUSSION 
3.1. Isotropic luminosity of masers 
Above, we found that water maser detection rate in- 
creases with [O III] A5007 luminosity (£[oiil]A5007)) ve- 
locity dispersion (cr), and optical luminosity (Mb). We 
speculate that these correlations may be a consequence 
of an underlying correlation of these parameters with wa- 
ter maser luminosity. The true luminosity of masers 
is difficult to measure be cause the ma ser emission is 
likely to be beamed (e.g., Elitzur 1992) and estimates 
of the beaming angle require a detailed model of the 
maser, which can only be inferred in the cases with well- 
understood geomet ries from VLBI observations (e.g., 
iMiyoshi et al.l 119951 ) . In place of true luminosity, we 
adopt apparent luminosity, which is based on the premise 
of "isotropic" emission of radiation. The isotropic lu- 
minosity can be computed readily from the flux densi- 
ties observed in spectra. On the other hand, the flux 
density can be va r iable on time scale s of months (e.g., 
Bragg et al.l 120001: IBraatz et al.l 120031: iHerrnstein et alJ 
20051: iCastangia et al.l 12008^ which introduces another 



uncertainty. Nonetheless, analysis using the stand-in of 
isotropic luminosity provides an opportunity to investi- 
gate whether the speculated correlations exist. From the 
literature, we have collected isotropic luminosities for 66 
masers. 

Maser surveys are flux-limited and the lower maser de- 
tection rate in galaxies with lower £[oin]A5007,obsj smaller 
a, and fainter Mb could be because most masers in these 
galaxies are too faint to be detected. We investigate the 
flux limit in Figure [9l where we plot isotropic luminosi- 
ties as a function of redshift. Although these detections 
are from a variety of surveys that have different sensitiv- 
ities, they appear to be consistent with 0.1 Jy km s _1 as 
the effective limit. This limit is consistent with a plau- 
sible detection threshold of 10 mJy (3er) in a 1 km s" 1 
channel and blends of Doppler components on the order 
of 10 km s~\ as seen in the spectra (e.g., with GBT, 
IBraatz et al1l200ll ). 

To investigate the correlations between maser luminos- 
ity and optical properties of the host galaxies, we comple- 
ment the sample with Mb, cr, i[oni]A5007,obsj and Ha/H/3 
(thus i[oni]A5007,cor) from the literature. In addition to 
those with SDSS photometry or spectroscopy, we have 
compiled L[ O iii]A5007,cor (Ha/H/3), £[om]A5007,obs, cr, and 
Mb for 36, 40, 43, and 54 masers with measured isotropic 
luminosities. We present these data in Table [TJ The 
combined sample is inhomogeneous in terms of selection, 
technique, and instrument parameters, and this (in ad- 
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Fig. 9. — Isotropic maser luminosity as a function of redshift 
for 66 extragalactic nuclear masers. The dotted line shows the flux 
limit 1.0 Jy km s , and the dashed line shows the flux limit 0.1 
Jy km s -1 . 



dition to the use of isotropic luminosity as a stand-in - 
noted earlier) may be expected to increase the scatter in 
any correlations. 



In Figure [101 



we 



log 10 £[OIH]A5007,obs 



log 



plot log 10 L 



10 ^[OIII]A5007,cor; 



10 1 



isotropic luminosity log 10 ^H20 



and Mb against the 
L-R20- Although the scatter 
in all relations is relatively large, the variables are 
clearly correlated. Assuming a uniform error of 0.5 dex 
for log 10 LH2Oj we perform least-squares fits with the 
following linear relations: log 10 £h20 = a + b (x — xq), 
where we choose xq to be the medians, 41.39 
(dex), 40.43 (dex), 2.16 (dex), and -20.58 (mag) 

for X = log 10 L[oiII]A5007,cor7 l°gl0 ^[OIII] A5007,obs, 

log 10 cr, and Mb, respectively, to minimize the cor- 
relation between a and b. We show the results with 
the dashed lines. The intercepts and slopes (a, b) for 
lo 

§10 ^[OIII]A5007,cor> log 10 -^[OIII] A5007,obsi l°gl0 °~ > an d 

M B are (1.65±0.13, 0.31±0.10), (1.66±0.12, 0.22±0.11), 
(1.64±0.11,2.67±0.64), and (1.75 ± 0.11, -0.21 ± 0.13), 
respectively. We can draw a robust conclusion that 
masers are stronger in hosts with higher [O III] A5007 lu- 
minosity, larger a, and higher optical luminosity, as 
implied from the detection efficiency results. 

We would like to emphasize that the best-fit slopes and 
intercepts here should be taken with caution. First, as 
mentioned above, the isotropic luminosity is a poor indi- 
cator of maser strength and the sample of optical prop- 
erties is not homogeneous. Second, stronger masers are 
more easily detected, and at the faint end there should be 
a detection bias favoring stronger masers scattered above 
the real correlations; the real slopes therefore should be 
steeper than the best-fit ones. Finally, the sample size 
is small and a few outliers could significantly affect the 
fitting. However, the conclusion is robust that maser 
strength increases with [O III] A5007 luminosity, cr, and 
optical luminosity. 

In Figure [H we showed that maser detection rate is 



Hosts of Water Masers 



7 



I 








m 


■ 




■ 


■ < J 




✓ 


■ 






■J 

- 












[■■■' > 




.-'/■ 










■ 


/ ■ ■ 






■ ■ 


■ / Fit: Liuo'lfoiii] • 


Fit: L Hao «L[o ! ; I ,] ■ 


■ / Fit: L H20 «(7 2 - 7 . 


Fit: L H20 «L°' 5 ■ 


' t: Lh30 k L[o°ii] 


Tf . T x l 1.0 

. 1 ... 1 ... 1 ... 1 ... 1 . 




■f: Lh30«Lb° 


38 40 43 44 


36 38 40 43 44 


1.4 1.6 1.8 3.0 3.3 3.4 3.6 - 


16 -18 -30 -33 -34 


lo gl0 L [01H] 15007. cor ( el "g S ~') 


lo gio L [om] A5007. ob S ( er g s ~') 


logio " ( km s ~') 


m b ( m ag) 



Fig. 10. — Relations between isotropic maser luminosity (I/H20) an d ^[Olll]A5007,cor> ^[Olll]A5007,obs> a i an< l Mb of the host galaxies. 
The dashed lines are the linear least-squares fits assuming a uniform error of 0.5 dex in log 10 Lh20- The dotted lines show the linear 
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O 



M 
O 




5 10 15 20 25 
H„/H« 



Fig. 11. — Relation between isotropic maser luminosity (Lh2o) 
and Ha/H/3 of the host galaxies. 

higher at higher Ha/H/3. As for the other parameters, 
this could result from an underlying relation between 
maser luminosity and extinction. We investigate such a 
relation in FigurcQTJ A linear least-squares fit gives (a, b) 
= (1.54±0.15, 0.03±0.03) for x = 6.31. We do not find 
a convincing relation (with only ~ la) between maser 
luminosity and Balmer decrement, indicating that there 
might not be a direct linear relation between log 10 Lh2Q 



and Balmer decrement. 
Section 13.31 



3.2. 



We discuss this result more in 



Correlations of maser emission with central black 
hole mass and AGN activity? 

It has long been proposed that water maser emission 
is closely related with the central black hole. Assum- 
ing that a thin viscous accretion d isk is obliquely illu- 
minat ed by a central X-ray source, iNeufeld fc Malonevl 
(1995) find that the critical outer radius R cr at which 
the disk becomes atomic and the maser emission ceases 



is the 



follows i? cr oc Lj°i 4 3 m a81 M^ 2 , where L 2 -io 
2 — 10 keV X-ray luminosity of the host galaxy, fa is 
the mass accretion rate, and Mbh is the central black 
hole mass. If we assume that the X-ray luminosity 
is proportional to AGN bolometric luminosity (Lagn), 
and that Laon is proportional to the accretion rate 
(e.g., [Frank, K ing, fc Rainel2 002). we consequently reach 
i? cr oc LS^bh 2 - Considering that in the geometrical 
maser models (e.g. iMivoshi et ailll995|) . maser spots do 
not cover the whole nuclear disk but rather lie on sev- 
eral radial arms, we assume that total mase r luminosity 
£h20 ^cr (but see lKondratko et aT]|2006bl who assume 
LH20 oc R 2 ct ). We then expect L H 20 oc £^n m bh 2 - If 
we further assume that AGNs radiate with a roughly 
fixed Eddington ratio r\ as a function of black hole mass, 
so that Lagn — V^Edd oc Mbh, we eventually reach 
Lh20 oc .Lagn oc M B h • 

It is well-known that there exists a tight corre- 
lation between the central black hole mass (Mbh) 
and the velocity d isp ersion of the black ho le host 
(jFerrarese fc MerrittJ l2000t iGebhardt et all l2000t 
iTremaine etHI [20021; IGultekin et alj|2009h : M BH oc cr 4 . 
The optical luminosity also correlates with the black 
hole mass but with a lar ger scatter ([Magorria n et al.l 
[19981: IGultekin et all [2009h : M BH oc L(Optical). For 
AGN strength, the [O III] A500 7 luminosity is a rea - 
sonably reliable indicator (e.g., iHeckman et al.l [2005) . 
In Figure [TQ1 assuming Lh20 oc Mbh oc Lb oc cr 4 and 
L mo oc £agn oc L[oiii]A5007, we fit the intercepts (at 
0) and obtain -39.6 ± 0.2, -38.6 ± 0.2, -7.0 ± 0.1, and 



-6.5 ± 0.1, for log 10 iv[ O III]A5007,cor, 



iO Sl0 ^[OIII]A5007,obs, 

log 10 cr, and Mb, respectively. We show these fits with 
dotted lines. If one takes our best-fit correlations in 
Section 3.1 at face value, it implies a weaker dependence 
of maser luminosity on black hole mass and AGN 
activity than the simplified theory. 

The analysis above is simplified. For example, the Ed- 
dington ratio rj is not necessarily independent of black 
hole mass for masers, even if it is so for AGNs. If masers 
are powered by AGN luminosity, then a high enough Ed- 
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dington ratio is required to pump up a strong maser. 
This hypothesis can explain what we found in Section 
12.2.21 that [O III] A5007 luminosity is more important a 
factor than Mb and a when determining detection effi- 
ciency. It is also interesting that the slopes above are all 
steeper than the best-fit ones. As we discussed in Section 
13.11 the best-fit slopes could be flatter than real ones due 
to the detection bias. Considering the large uncertainties 
caused by the variability and the assumption of isotropy 
in maser luminosity and the small sample size, we con- 
clude that our fitting results are in reasonable agreement 
with the simplified theory and that maser strength is in- 
deed correlated with the central black hole mass and the 
AGN activity of the host galaxies. 

Previous studies have also found correlations be- 
tween the isotropic luminosity and va r ious pr operties of 
the host galaxies. iKondratko et al. (2006b) find that 
£h20 oc^ £2-101 where L2-10 is the 2 — 10 keV X-ray lu- 



minosity of the hos t galaxy. iZhang et al.l (20101 see also 
Henkel et all 120051: iCastangia et al.l 120081: iSurcis et al.l 



2009) find there appears to be a correlation between 



Lfir and Lfj20i where Lfir is the total far- infrared 
(FIR) luminosity of the host galaxy. However, they 
do not claim to find a significant correlation between 
extinction-corrected £[oiniA5007,cor and Lh20: likely due 
to the smaller sample size in their study. 

3.3. Masers favor high- extinction systems: a 
geometrical effect? 

In Figure H] we showed that maser detection rate 
is higher at h i gher Ha/1 10, thus higher e xtinction. 
IGreenhill etHI (f2008h and IZhang et alfpOloh also find 
that there is a high incidence of Compton-thick systems 
among AGN masers. If masers are only located on the 
nuclear disk surrounding the central black hole, then 
this preference of higher-extinction systems could be ex- 
plained as a geometrical effect. 

The specific intensity of maser emission scales with the 
cube of the pump ing gain l ength (I) along the path (e.g., 
IStrelnitskiil[l98l OE005): Lmo °c I ■ If masers orig- 
inate in gas clouds in the circumnuclear disk, we can 
only observe masers when the disk is inclined so that 
the gain length along the line-of-sight is long enough. 
According to the propos ed unified AGN models (e.g., 
lAntonucci Sz. Millerl 119851 ). when the circumnuclear disk 
is inclined, the optical extinction due to the dusty torus 
and the X-ray obscuring column (iVjj) are correspond- 
ingly higher. Naturally, masers are more often found to 
be in high-extinction systems. Under this theory, we 
expect total isotropic maser luminosity on average to 
be higher in higher-extinction systems, since the aver- 
age gain length should be longer. However, we did not 
find a convincing linear relation between log 10 Lh20 and 
Ha/H/3 in Figure [TT] This could be because Ha/H/3 is 
not a direct proxy of gain length, or the intrinsic scatter 
of the relation is too big to detect in the small sam- 
ple we curr ently have. Meanwhi le, some masers, e.g., 
NGC 1052 dClaussen et alJ H99l) and NGC 262 (Mrk 
348, IPeck et al.ll2003D . are likely to be associated with 
jets but not disks. These systems could have greatly in- 
creased the scatter of the relation. It is also possible 
that the inclination of the host galaxy can affect the ob- 
served Balmer decrement, and thus increase the scatter. 



On the other hand, if this theory is correct, then many 
extragalactic nuclear masers should be disk systems, but 
more than half of current maser detections lack evidence 
of association with disks (or jets). To confirm this the- 
ory requires further investigations with larger samples 
and follow-up VLBI observations of maser detections. 

4. SUGGESTIONS ON SURVEY STRATEGIES 

Our results suggest that if we can improve the obser- 
vational sensitivity, we should be able to detect more 
nuclear masers. Given the current sensitivity, however, 
we can still improve the detection efficiency. Since there 
exist large optical spectroscopic surveys, such as SDSS, 
2dFGRS, and 6dFGS, we suggest that maser surveys pri- 
marily interested in efficiency should select AGN tar- 
gets from these surveys and rank them by extinction- 
corrected [O III] A5007 flux. As an example, we illustrate 
this strategy in Figure [TOl 

In Figure [H we plot log 10 L[oniU5007 against red- 
shift for all the 33 detections and 1030 non-detections 
with SDSS spectroscopy. Note in Figure [5] we showed 
that 0.1 Jy km s" 1 represents the effective flux limit 
in maser surveys. To relate this limit to £[oiii]A5007j 
we use the dotted lines in Figure [TO] log 10 Lh~20 = 
-39.6 + log 10 L[oni]A5007,coi and l°gio L H20 = -38.6 + 
lo gio £[oni]A5007,obs, where L H 20 is in L© and L[oiii]A5007 
is in erg s _1 . This conversion translates the limit 0.1 
Jy km s" 1 in maser flux into 7.6 x 10 -14 erg s _1 cm~ 2 
in extinction-corrected [O III] A5007 flux, and 7.6 x 
1Q-15 erg s -i cm -2 in b S erved [O III] A5007 flux. 
These limits are the solid lines in Figure [TO] To take 
into account the scatter in the £H2O--k[oni]A5007 re- 
lation, we shift the limits downward by —1.0 dex to 
7.6 x 10~ 15 erg s _1 cm~ 2 and 7.6 x 10~ 16 erg s -1 cm~ 2 
and show them with the dashed lines. Among the 33 
detections, 30 (32) have i[ O ni]A5007,cor (£[om]A5007,obs) 
brighter than the dashed line. For £[om]A5007,cor 
(-^[Oin]A5007,obs)) the detection rates above and below 
the dashed 'line are 6.9% ± 1.2% (4.3% ± 0.8%) and 
0.5% ± 0.3% (0.3% ± 0.3%); the detection rates above 
and below the solid line are 16.0% ± 4.1% (9.2% ± 2.2%) 
and 2.0 ± 0.5% (1.9% ± 0.5%). 

Therefore, if we rank AGN targets by extinction- 
corrected [O III] A5007 flux, we can significantly im- 
prove the detection efficiency, from an overall ~ 3% to 
- 16% for the strongest [O III] A5007 emitters. Fur- 
thermore, given the current sensitivity, observing AGNs 
with extinction corrected [O III] A5007 flux lower than 
the dashed line yields a very low detection rate of < 1%. 
The contrast among these detection rates suggests that 
ranking source lists according to extinction-corrected 
[O III] A5007 is effective in maximizing detection effi- 
ciency. If extinction correction (i.e., Ha/H/3) is not avail- 
able, we suggest that maser surveys rank AGN targets 
with the observed [O III] A5007 flux. If the observed 
[O III] A5007 flux is not available either, we suggest that 
maser surveys rank AGN targets by velocity dispersion, 
or by optical luminosity. This strategy should give a 
higher detection efficiency than a blind survey. 

We are grateful to Jim Braatz for allowing us to 
use the compilation of the galaxy sample surveyed for 
water maser emission before publication. We thank 
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Fig. 12. — [O III] A5007 luminosity as a function of redshift. We plot all detections and non-detections in the SDSS spectroscopic sample. 
The solid lines show the assumption of maser flux limit 0.1 Jy km s _1 and Lh20 x ^[OIII]A5007 (t ne dotted lines in Figure flOt . The dashed 
lines are the flux limit shifted downward by —1.0 dex, to take into account the scatter in the I/H20 _ ^'[0III] A5007 relation. The percentages at 
the end of arrows show the detection rates above the corresponding lines and the two numbers show the number of detections and galaxies 
surveyed for maser emission (detections plus non-detections). 
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TABLE 1 

Physical parameters of extragalactic water maser sources 



Source 


V 1 

v sys 


LH20 2 


Rcf(L H2 o) 3 


Froim 4 




Ref(F [on i]) 6 


a 7 


Rcf(<r) 8 


Ms 9 


Ref(M B ) 10 


NGC 23 (MrK 545) 


4566 


2.3 


BG08 


1200 


5.2 


DD88 






-21.50 


DeV91 


NGC 17 (Mrk 938) + 


5881 




















2MASX J00114518-0054303 


14384 
















-19.58 


SDSS 


NGC 235A 


6664 


2.0 


Kondratko06 








198 


Wegncr03 


-21.02 


DeV91 


NGC 262 (Mrk 348) 


4507 


2.6 


Henkel05 


22710 


6.02 


Bassani99 


185 


McElroy95 


-20.54 


DcV91 


NGC 291 


5705 






6431 


4.97 


SDSS 


109 


SDSS 


-20.14 


SDSS 


ESO 013-G012 


5047 




















NGC 449 (Mrk 1) + 


4780 


1.7 


Henkel05 


60000 




Whittle92 


115 


NW95 


-19.53 


DeV91 


NGC 520 


2281 


0.1 


Castangia08 


412.4 


4.13 


Ho97 


41 


Ho09 


-20.55 


DeV91 


2MASX J01260163-0417564 


5639 




















NGC 591 (MrK 1157)+ 


4547 


1.4 


Henkel05 


23000 




Whittle92 


95 


NW95 


-20.50 


DeV91 


NGC 613 


1481 


1.2 


Kondratko06 








125 


McElroy95 


-21.08 


DeV91 






1.2 


Castangia08 














IC 0184 


5382 


1.4 


Kondratko06 












-20.00 


DcV91 


zivi/\ojV j uzi4uoy i-uuioo 1 1 


1 1 one; 

1 1ZU0 
















1 fi e;o 

-lo.oy 


olJoo 


ivirK luzy 


yu (D 








"' , 












NGC 1052 


1510 


2.1 


Hcnkcl05 


• > » 1 / WW) 

33068 


2.35 


Ho97 


187 


Wegner03 


-19.85 


SDSS 


NGC 1068 (M 77)t 


1137 


2.2 


Henkel05 


3497963 


5.29 


Ho97 


162 


Ho09 


-21.68 


DeV91 




1.7 


KGM06 
















NGC 1106 


4337 


0.9 


BG08 








146 


Wegner03 


-21.17 


DeV91 


2MASX J02532956-0014052 


8622 






1757 


5.0 


SDSS 


96 


SDSS 


-18.18 


SDSS 


Mrk 1066 


3605 


1.5 


Hcnkel05 


24000 


8.51 


Whittle92 


105 


NW95 


-20.60 


DcV91 


NGC 1194+ 


4076 






2340 


5.85 


SDSS 


144 


SDSS 


-20.06 


SDSS 


NGC 1320 (MrK 607) 1 " 


2663 




















2MASX J03364614-0750236 


11719 




















NGC 1386+ 


868 


2.1 


Henkel05 


99244 


5.7 


Bassani99 


187 


McElroy95 


-18.52 


DeV91 


2MASX J03381036+01 14178 


11926 


















(IRAS03355+0104) 






















4C +05.19 


790000 




















(2MASX J04143774+0534423) 
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TABLE 1 

Physical parameters of extragalactic water maser sources 
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TABLE 1 

Physical parameters of extragalactic water maser sources 



^ Heliocentric systemic velocity, in km s _1 . 

Total maser luminosity assuming isotropic emission of radiation, defined as Lh20 — 0.023 X f S(V)dV X D 2 , where Lh20 is in Lq, S(V) is the 
flux (in Jy) at velocity V (in km s" 1 ), and D is luminosity distance (in Mpc). These values were calculated based on Hq — 75 km s~ 1 and We 
convert values into those with Ho — 70 km s -1 in the analysis. 

3 References - Henkel05: IHenkel et~aL1 1 )20051) ; Kondratko05: IKondratko et all d2005l) ; Kondratko06: IKondratko et al.| d2006a|) ; KGM06: 
iKondratko et all d2006bl) : Zhang06: iZhang et all d2006D : BG08: iBraatz Gugliuccil d20Qgft : Castangia08: ICastangia et al.l d20Qgft : Surcis09: 
ISurcis et all d2009l) . 

4 Observed flux of [O in] A5007, in 10~ 17 erg s" 1 cm" 2 . 
Ratio of observed flux of Ha to that of H/3. 

6 References - DD8S: IDahari & De Robertlsl JT988T) ; Whittle92: IWhittlel dl992ft ; Polletta96: IPolletta et~aT1 dl996D ; Ho97: I Ho et al.| fl997l) ; 
Bassani99: Bassani et al. (1999). In the compilation, we first prefer the homogeneous samples from SDSS and Ho97. We then choose the most 
recent values for the rest of the sample. 
Velocity Dispersion <r, in km s — . 

8 References - McElrov95: iMcElrovl fl995h : NW95: [Nelson & Whittle! fl995h : Qliva95: fOliva et al.l fl995h : Qliva99: fOliva et al.Nl999H : Wegner03: 

IWegner et"aT1 d2003l) : Gerssen04: IGerssen et~aT1 J200i) ; Gu06: IGu et all d2006fl ; GH06: IGreene & Hoi J2CJ06I) ; Ho09: IHo et all fl2009l) . 

9 

B band absolute magnitude. 

References — DcV91: dc Vaucoulcurs et al. (1991). Magnitudes are corrected for foreground Galatic extinction ijSchlcgcl et al1 |1998T) . 
a The merging system Arp 299 consists of two galaxies, NGC 3690 and IC 694. We assume the maser is associated with NGC 3690. 
' Assumed to be 3 in the analysis. 
' Likely disk systems. 



